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Note on certain Coefficients appearing in the Algebraical 
Development of the Perturbative Function. By B. T. A. Innes. 

{Originally written in January 1909 {but lost in post), and rewritten 

1 8 th. March 1909.) 

I. 

We find a short algebraical development of the Perturbative 
Function in terms of the mean anomalies in Laplace’s Mecanique 
Celeste. This was extended by Burckhardt and de Pont^eoulant. 
Later, Le Yerrier gave a development to the 7th. powers of the 
eccentricities and 6th. powers of the inclinations (both inclusive), 
and this was extended to the 8th. powers by Bouquet. Le Yerrier 
gave extended tables and numerous precepts with his develop¬ 
ment ( Les Annales , i. and x.), and in consequence it has been 
widely used. 

A more methodical and general development was given by 
Newcomb (Ast. Papers , v., parts 1 and 4, 1895). Earlier, New¬ 
comb had published a development in terms of the eccentric 
anomalies (Ast. Papers , iii.) which he used in deriving the 
periodic perturbations of the four inner planets, but it would seem 
that the further adoption of this development has been tacitly 
condemned by its illustrious author—it appears to be necessary to 
state this explicitly, as in a recent didactive treatise on this subject 
(I refer to Mons. H. Poincare’s Legons de Mec. Cel., ii., part 1, 
1907) the author’s references to Newcomb omit his later and more 
generally useful papers. 

The earlier developments (including Le Yerrier’s) make use of 
the Laplace coefficients 

b^\ b^\ b^, etc. etc., 

and of their derivatives of the first and higher orders. An in¬ 
spection of Le Yerrier’s formulae will show that the coefficients 
which actually appear are in a slightly different form, viz. 

bf, ab^\ a?b etc. etc., 

and that the derivatives of these are there implicitly. 

Finding the numerical values of the various Laplace coefficients 
(especially if a is over J and accuracy is required) is a wearisome 
process, but when a is less than 075 the actual work required can 
be lessened by using tables by Bunkle and others. Well-known 
recurrence formulae exist between the various coefficients, but un¬ 
fortunately their use leads to accumulating errors.* Le Yerrier 

* This was pointed out by Le Yerrier in one of his earliest papers (Paris, 
1841) and is well known, yet we find Mons. H. Poincar 4 , in the work 
already cited, recommending their use. Of course the recurrence formulae are 
mathematically exact, but they are too inaccurate for astronomical use, 
paradoxical as it sounds ! 
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introduced transformations of power-series so as to facilitate their 
calculation, and gave extensive tables (Les Annales , ii.), but even 
Le Verrier’s figures will not give the 5th. place correctly if a is so 
large as 0-825 (Titan-Hyperion). 

Besides the recurrence formulae, Le Verrier gave {Les Annales , 
ii. pp. 108-9) several isolated relations between groups of the 
coefficients, such as 




da 


do? 


He found ten such relations, but gave no indication that a 
generalisation exists. The advantage of such relations is that they 
consist of simple additions of tabulated quantities multiplied by 
integers. The recurrence formulae previously mentioned involve 

awkward multiplications by fractional numbers and by — + a, 

a 

1 + no , 2 , etc. 


II. 

The object of this paper is to give a general formula which 
leads to a very simple recurrence formula, and embraces the 
generalisation mentioned above. This will be done for Newcomb’s 
development. I will use the symbol introduced recently by 
Cowell—it has the great advantage of avoiding the use of an affix 
in the usual place of a power-index. The coefficients will compare 
as follows :— 


Laplace. 

Newcomb. 

This paper. 

bf 


hi 

a bf 

r (i) 

C 3 

h 3 


Ji) 

C 5 

h 5 

etc. 

etc. 

etc. 


In place of the ordinary derivatives, Newcomb has introduced 
logarithmic derivatives, indicated by the operator D = a ^~. The 
general relation is 

(«|J=D(D-i)(D-2) .... (D -n+ 1), 

and it enables us to pass from one set of derivatives to the other 
if wanted. 

The recurrence formula used by Newcomb 
m,j)_ n+V+ no? „ (ilJ+1) _,_(» + 2 j)(zi +2 -n) 

W ~ ITpj C * + ” 

is free from some of the defects of that of Laplace. The quantity 
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c { n ij) is an intermediary between and c%\ which vras intro¬ 

duced by Newcomb to facilitate calculation.* It is, to a simple 
factor, a hypergeometric series. If we apply the Gaussian 
relation 

F/3+ - F = — vF'" 

7 

to this quantity we get 

(n - if 4* +1,i) = (2W +j - 1 )4-i +1) + 4-i +2) , 

a recurrence formula freed from all defects. But what we really 
require are the quantities c {i) , Dc (i) , D 2 c (i) , etc. If we apply the 
formulae given by Newcomb to find these quantities (Ast. Papers , 
v. P* 3 1 3 j etc.) the last recurrence formula becomes 

or using the newer notation 

(» - *fV m h, n = D“(d + i + ^ 3 )(d - i -^ 5 )^ 1 .. 

Particular cases are— 

Dm c (0) = Dw p 2+ D-o)4 1) 

D m 4 P = B m (D 2 + D - o) 4 0) = D m (D 2 + D - 2 ) 4 2 > 

D m 4 2) = D m (D 2 + D - 2 ) 4 j) = D W (D 2 + D - 6)cf 


(n - 2) 2 D m 4 i) = D m ( D + i + 


9 D m 4 0) = D m (D 2 H- D — o)4 1) 

9 I) m cU = D W (D 2 + D - 2)4 0) = D m (D 2 + D - o) 4 2 > 

2 5 D m 4 0) = D m (D 2 + D - 2 ) 4 11 
2 5 D-4D = D m (D 2 + D - 6 ) 4 0) = D W (D 2 + D - o)<f 

4 9 D w 4 0) = P m (D 2 + D - 6 ) 4 15 etc. etc. 

With these new formulae we may check some of the figures 
furnished by Newcomb (Ast. Papers, v. p. 343) for Venus-Earth. 


We have 

mf = 

66839 

(Newcomb) 

Check 

D 5 <f = 

646 9 i 



D 6 <f = 
-6D 4 cf = 

4182 
~ 2033*4 



Hence 

6683 9 *6 = 

:D 4 4 8 ». 


* It was also used in the same sense by W. Soheibner (Storungsf unction. 
Gotha, 1853) 
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If I> 4 4 4) is thus calculated, it is found =65583, whereas 
Newcomb has 65850, but an independent calculation has refound 
the former figure. 

Again, 

D 4 4) = 60000 Newcomb (apparently the 

- nearest round figure 

Check D 3 4 3) = 2974000 which is sufficiently 

D 2 „ = 118490 accurate). 

-42D „ = - 224952 

1/49 x 2867538 = 58521 =D4 4) . 


From the general formula many simple relations may be 
deduced, such as 

D ”‘(*0,3- & 2,s) = 2l)m *l,l 

3D"‘(A), 5 — ^2, 5) = 2 ^1, 3 

5 Dra ( 6 o, 1 “ K1) = 2D m K 5 etc. etc. 


The importance of the general formula in the development of 
the Perturbative Function is considerable. It not only simplifies 
the calculation of b i 3i b i5 , etc., but it can also be turned to use in 
computing the D m b { L coefficients, which will save the introduction 
of the quantities. 

We have Newcomb ( Ast. Papers , v. p. 310), and noting that 

D4 i) = ^f ) + 4 i * 1) 


i(D -i)b itl 


1 * 3-5 


2fc+ I 


2.4.6 


21 + 2 


a <+2 F(f, 1 +?;, 2 + 3 , a 2 ) 


which transforms into 


2.4.6 . . . 21+2 

wherein 


• • 2 <+I «‘^F(i, J + i, 2 + 3 , a 2 ) 


/ 3 2 = —„• 

I — a/- 


Tisserand (Mec, Cel., i. p. 274) 




J 


shows that 

sin 2 * 0 dQ 
1 - a 2 sin 2 0 


and we also have 


F(i, i + i 9 2+t, a 2 ) = N— 


sin 2 * 0 - sin 2 * +2 0 
Ji - a 2 sin 2 0 


■d0, 


wherein N is the ratio of two whole numbers. 
Hence 

(D - i )\! = (23 + 

and x can be found from b it : and b i+1> 
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Later, it will be seen that we should compute 



and this must be done with precision to avoid the introduction of 
error. 

To find the higher derivatives we can use the differential 
equation satisfied by the coefficients (Tisserand, Mec. Cel i. p. 
281) which in the new notation is as follows— 


(D 2 - i 2 )b it1 = / 3 2 ( 2 D + i ,, 


and noting that 

D/3 2 = 2 /5 2 (i + / 3 2 ) 

we have 

D(D 2 - i 2 )b it ! = /i 2 [4D 2 + 5D + (2 - 2 i 2 )]b it x 

D 2 (D 2 - i 2 )\ x = / 3 2 [ 6 I) 3 + 13D 2 + (12 - 4* 2 )D + (4 - 4 i 2 )]K i, 

etc. etc. 

To obtain either D m+ 2 b 0>1 or D m+2 6 L1 we must use these 
formulae; for other values of i there is a better process. Between 
D^o, j and T) m b lt 1 we have the following simple relations — 

D6 0i 1 = aDfr L j 

D 2 6 0il -=aD(D + i)6 la 


D m b 0i x = aD(D + 1 ) m ~ 1 b h ] . 

To obtain D m+ 2 b itl for other values of i we can write the 
formula— 


D m b 


i, n 



D 2 + D- 



w -3 


b 


t-1, n-2 


first with i positive and then with i negative. Equating the two 
values so found and confining our attention to the case of n = i 
yields 

D-(D 2 + D-2)6 2>1 = D™(D 2 + D-o )6 0§1 

D m (D 2 + D - 6 )b Bt j = D W (D 2 + D-^2)& L1 



etc. 


etc. 


in actual work we have 




D” + \, 

= D”‘(D 2 + IJ 

1 

O 

© 

- D“(D - 

■Oki 

D” +2 & 3 , j 

= D m (D 2 + D 

~ 2 )\, 1 

- D"*(D - 

~ 6 )b s , 1 

D" +2 & 4i x 

= D m (D 2 + D 

“ ^)^2, 1 

- D m (D - 

- I 2 ) 6 4 ,I 


etc. 


etc. 
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Thus all the quantities except the b itl can be derived quite simply 
once the latter are found, no large numbers nor awkward factors 
are involved, and there is no accumulation of error except that 
which is inevitable in any extended calculation, and which arises 
from the uncertainty of the last decimal place. 

The formulas developed above are useful in computing the 
action of the planets on the Moon, and evade a difficulty to which 
Brown has called attention (Action of Planets , p. 34). 

It will be apparent that the introduction by Newcomb of the 

operator 13 instead of a~ is not only an analytical improvement, but 

da 

one of great practical importance as well. 

I am now finishing a numerical application to the case of Titan- 
Hyperion, for which a exceeds 0*855, in which the whole process, 
with checks and alternative methods, is fully illustrated. 

Johannesburg , S.A. 


Numerical Example. (Innes.) 

Calculation of b 0i v 

This is the most important constant amongst the coefficients 
1 n , as its value enters into all of them. 

Earlier investigators computed its value by means of the hyper- 
geometrical series F(£, 1, a 2 ), which is not very converging, 

especially for a>0*707. Le Yerrier gave a transformation of this 
series (Annates, ii. p. [5]), which was used by Runkle in his well- 
known tables, but the transformation is still unsuited for a rapid 
and exact calculation, and in fact Runkle’s 7th. decimal place is 
frequently in errqr; thus for a = 0*75 Runkle gives 

log J 01 = 0-3861679, 
whereas the true value is =0*3861685. 

Newcomb, in his planetary theories (Astr. Papers , iii. p. 69), 
recommends a formula suggested by the use of Gauss’s arithmetico- 
geometrical mean, viz. 

log b 0l — N - log cos - £ log cos 0 (a = sin 0 ) 

wherein N is very nearly equal to log 2. A small table of N is 
given by Newcomb, but it has evidently been computed by means 
of Runkle’s table, as it suffers from the same errors. The value of 
N is probably equal to 

log 2 + 2 log(i - 2(f + 2g 16 - etc.) 
wherein q is the quantity known as Jacobi’s Nome. 
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I have preferred to calculate the value of b Q> j by means of the 
following formulae, which have been adapted from Verhulst 
(Fonctions Elliptiques , Bruxelles, 1841, p. 256):— 

8 

K 1 = ( I+ < y-)2 ( I + 2 g i +2q w + etc.) 2 
or 



1 + 2£ 16 + 2g 64 + etc.) 2 . 


It seems inadvisable to replace cq by cos 0 . 

In the first of these two formulae, the terms in q amount to less 
than one unit in the 8th. place of decimals if a does not exceed 
0*35; in the second formula the terms in q will only amount to 
one unit in the 8th. place of decimals if a is as large as 0*998, a 
value much larger than is ever met with. In the case of Titan- 
Hyperion (a = 0*825), the use of the second formula, neglecting the 
g terms, introduces an error of 7 in the 19 th. decimal place of the 
logarithm of b 0t v 

If we write 


lp g K i = log 8 - 2 log (1 + JaJ + Nj 

the quantity N 1 may be taken from the table on p. 649. 

The following values of log q for different values of a will afford 
a rough estimate of the q terms in b 0t 1 :— 


a = 0*05 
0*15 
0*25 


log q = 6*20 

« = o*35 

l°g ? = 7 ' 9 I 

£, = 0*65 

7 * i 5 

o *45 

8-i 5 

075 

7*6i 

o *55 

8'35 

0*85 


{N.B .—10 has been added to log q .) 


log q = 2 >'53 

871 

8-90 


As a numerical example we compute log & 0(1 for a = 0*75 by 
the approximate formula 

log b o, 1 = log 8 - 2 log j~ x + ^ a i + sf il t ^ a l . 



log cq = 9*82048903 


1 + s/cg 

2 


0*90664414 



(i + a i) Vcq 
2 


0*90661866 


log 8 = 0*90308999 
2 log 1*81326280 = 0*51692150 

log V i = 0 ‘3 86i68 49 


It will be noticed that although q has been introduced into the 
formulae, it is in reality not used (owing to the extraordinarily close 
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approximation to the true value given by the second formula) ; thus 
it does not seem necessary to repeat the equations by which q may 
be found from a; they will be found on p. 494 of the M.N., 
vol. lxii., where, however, \ is written for cq. 

The next step is to calculate b it ! from i— 1 to such a value of i 
as experience may indicate. Having found b 0t l5 there is no better 
way of finding the other values than that first shown by Hansen, 
and more fully and simply set out by Newcomb ( Ast. Papers , iii. 
p. 64). The process here followed only deviates in unessentials 
from that of Newcomb. I prefer to write 

= fro, 1 Pi P2 Pb 

a 

because the quantities p change slowly for considerable changes 
in a. At first sight the computation of these quantities appears 
complicated, but on closer acquaintance it is found most simple. 
Once it is started, the only table required is Zech’s subtraction- 
table, and in general only one opening of that. The one dis¬ 
advantage is that Zech’s table has at most only 7 decimal places, 
and no more extended table is available; hence it is not possible 
to be certain of the last figure; in spite of this, the exactitude is 
quite sufficient in practice. 

The quantities b 0>1 and p v p g, . ... p t are computed so 
quickly and easily that tables giving them with a as argument 
seem to be superfluous. To exhibit, however, the march of these 
quantities, Mr. Frank Robbins has computed the short table 
on p. 650. Short as this table is, it will permit the p i to be 
taken out almost correctly by inspection. 

In the case of Titan-Hyperion, for which a — 0*8250634 (H. 
Struve), the following are the chief figures involved— 

log <x = 9*91648732 

1 -<*^=<^ = 0-31927039 log = 9*50415864 

i + \Ah = 1 *75 i6 9 i 54 l°g = o*2 4 34576 3 

Hence by the formulae on p. 639, 

fr 0)1 = 2-6074655 ^ = 0*41621857 

To calculate p n use is made of the well-known Gaussian continued- 
fraction. In this calculation the letters a, b, c, d . . . . and F are 
simple fractional factors; Zech is put against the figures for which 
Zech’s subtraction-table is used. The calculation is quite similar 
to the example given by Newcomb Papers , iii. p. 71), but is 

more extended because of the large value of a. 


co-log 

Zech 

/ 9-83 

O 2 9 '83 

- Q ’34 

0*265 
e 8*416 

9-833 

- 1-486 

0*0144 

d 9*8716 
9*8330 

- 0*2803 

0-32199 

c 8*09211 
9-83297 

- 175292 

0*007740 

b 9*928271 
9-832975 

- 0*231014 

0-3845441 
a 7-2773661 
9-8329746 

- 2-5051152 

0-0013594 

F 9 * 97979 66 

9-66 -1 8-514 - 

97197 - 

8*24708 - 

9768986 - 

7-4948848 - 

9-9811560=7?!! 
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With this value of p Y1 we proceed 


Pn = 

9*9811560 

a 2 /( i + a 2 ) 

9*6074768 

Constant 

0*0211893 


9*6098221 

co-log 

°' 39 OI 779 

Zech 

0’227IOI7 

Constant 

9-9777236 

1/(1 +a 2 ) 

9 ' 7745 021 

Pi o . 

9 ’ 9793 2 74 


and so on; the only tabular figure in each block being taken from 
Zech. 

The following values of p result— 


log p 1 = 

9-7560111 

P2 

9-8960884 

Ps 

9-9322424 

Pa 

9-949308° 

P5 

9 ' 959353 ^ 

Pe 

9*9660080 

Pr 

9 ' 97 ° 755 6 

Ps 

9 ' 9743 201 

P9 

9-9770987 

P10 

9-9793274 

Pn 

9*98l I560 


In possession of these numbers we proceed to find 

^1,1 = 1 Pi 

& 2t 1 = a b\, 1 P% etc. etc., 

then Db it 1 by 

(D-o)& 0 , 1= = 

( D - 1)^1.1 = 3 ^* 1 , i( 1 ~Pi) 

(D - 2)fi 2i 1 = 5 / 3 2 & 2 ,1(1 -p s ) 
etc. etc. 

We have now to find D' 2 +m b Qt j and D 2+,n & L l5 which do not fall 
into a more simple general formula for D 2 +m b i+2jl . We therefore 
use 
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D°(D 2 - i 2 )\! = /3 2 [2D + i\b it ± N.B .— b it x is omitted in the next lines. 

D(D 2 -^ 2 ) = /3 2 [ 4 D 2 + 5 D + 2(1 - a 2 )] 

D 2 (D 2 - i 2 ) = /3 2 [6 D 3 + 13D 2 + (12 - 4 * 2 )D + 4(1- i 2 )] 

D 3 (D 2 - i 2 ) = /3 2 [8D 4 + 25D 3 + (38 - 6* 2 )D 2 + (28 - 12t 2 )D + 8(1 - i 2 )] 

D 4 (D 2 - i 2 ) = /3 2 [ioD 5 + 4 iD 4 + (88 - 8* 2 )D 3 + (104 - 2 4 * 2 )D 2 

+ (64 - 32i 2 )D 4- 16(1 - i 2 )] 

D 5 (.D 2 -i 2 ) = yS 2 [i 2D 6 + 61D 6 + (170— icn 2 )D 4 + (280 - 4 o2 2 )D 3 

+ (272 - 8o& 2 )D 2 + (i 44 - 8oi 2 )D 4- 32(1 - i 2 )] 
D 6 (D 2 - i 2 ) = /3 2 [i 4 D^ 4- 85D 6 4- (292 - 1 2 * 2 )D 5 + (620 - 6o^ 2 )D 4 

4- (832 - i6oz 2 )D 3 + (688 - 24.cn 2 )D 2 4- (320 - i922 2 )D 4- 64.(1 - i 2 )\ 

Having found D m & 0>1 and D m b lt 1 by these formulae we may 
check the values by means of the following equations:— 

D6 0> j = aDb lt x 
D(D - 1 )b 0f j — aD 2 6 L 1 


J)(D-i) m b 0tl = aD m+ 1 b hl . 


When satisfied as to the accuracy of the values obtained for 
these coefficients, we proceed with 


V m+ \. 1 = D TO [(D 2 + D - o)J 0 , j - (D - 2)b % J 
V m+ \ 1 = D“[(D* + D- 2 )b h , - (D - 6 )b 3 , J 


D m+ \ 1 = D“[(D 2 + D - (i- i)(» - # Hi - (D - i(i- i))b t , J. 


With the above sets of formulae the following table has been 
computed:— 


Titan-Hyperion (a = 0’8250634). 


i 

hi 

D&i,! 

&hi 

1 

D%, 1 

D5&i,l 

v*h 1 

0 

2*607466 

2-389585 

1574935 

170*9130 

2706*378 

56736-48 

1482263 

I 

1*226640 

2*896245 

16*19242 

171*8703 

2713*127 

56814-88 

1483538 

2 

0796696 

2*820306 

l6*9I202 

175*3909 

2735725 

57057-92 

1487413 

3 

0-562370 

2*611792 

I7*39780 

180 5431 

2776*454 

57491*06 

1494094 

4 

0*412874 

2‘3593 6 ° 

I7*54728 

186*1462 

2834*065 

58140*99 

1503924 

5 

0*310212 

2 *098800 

17*36659 

191*2088 

2904*347 

59020*83 

1517334 

6 

0*236676 

1 -847254 

l6*902l8 

195 0217 

2981*309 

60121 *47 

1534702 

7 

0*182556 

1-613119 

l6*2I328 

197*1492 

3058*366 

61410*81 

1556265 

8 

0*141972 

I 400109 

15‘35940 

197 3859 

3129*180 

62836*30 

158200*6 

9 

10 

11 

0*111120 

0*087419 

0*069063 

I *209330 
I*04039I 

H‘39452 

1957050 

3188*272 

6433i‘3i 

1611631 


12 (o’05478 est.) 

13 (0*04360 est.) 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at University of North Dakota on May 25, 2015 




1909MNRAS..69..6331 


June 1909. Development of the Perturbative Function. 643 
For values of D m b it7l (n other than 1) the formulae 

n+2 =D»*[D2+D - (i + + n -^fj\h » 



will be used, firstly with n — 1, then %= 3, 5, 7, and so on. 

The following empirical formula, which is only valid for large 
values of a, is curious— 


D<\ x = D«6 0i , + i[D« - D« + D*%, r 


In the case of Titan-Hyperion this shows no deviation to 7 places 
of figures. 

It will be remarked that from first to last no term of an 
infinite series is computed; the nearest approach to an infinite 
series is the continued fraction for_p u . 

It is well to have checks of an entirely independent nature on 
the values of the coefficients. Suitable checks will be found in 
the following formulae, where indeed use is made of infinite series, 
but these series are so convergent that even with a so large as 
0*825 . . . three or four terms are sufficient to ensure accuracy 
in the eighth significant figure. 

These formulae have been derived by Gaussian transformations 
of the ordinary hypergometric series. 

Put 

a = sin 0, 


and to abbreviate 

Q 

sin 0 tan 2 0 = st and tan 4 — = x ; 

2 

then 

* 0 , 1 = 2 ^j F (h b *• *) 
cos 2 — 

2 


a x — cos 0, — — tan 0, 


9 0 I + cu 

cos J — =-h 


(D +1)60,! = 


(D+ i)(D + 3 )6 0 j = 


,0 

COS'*— 

I.*) 

tt l 

<9 


COS 6 — 

1.3.2f, -I, I. as) 

a l 
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6 


cos 


10 _ 


(D+i)(D + 3)(D + 5)6 0i1 = I-3-S-2 „ F(-|, i.x) 

0 


cos 


14 _ 


(D + i)(D + 3) ... (D + 7 )b Qt 1 = I.3.5.7. 2-— t-F( - f, - h ”1. x) 

U-1 


(D- i)J u =- 3 .28 t 1 i, 2- *) 


2.4 0 _ 

^ COS 2 — 


e_ 

2 


cos - 4 

(D + 2)(D- i)6 lfl = |-|.2st—-i, 2. x ) 


(D + 2)(D + 4)(D - i)&j, 1 = 2 st — ^F( - 4, 2.x) 


cos' 


<9 


6 __ 


(D - 2)(D - 4 )6 2 , x = 1.3. “g 2 (st)* ^F( - 1 , i, 3 - *) 

2 

6> 


+ 3 )(^ “ 2 )(-^ " 4 )^ 2 , i — T, 3 « 2 .4.6.8 2 ( S ^) 2 a \ ^ h 3 - x ) 

J 

cos 6 - 

( D + 3 )( D + S)P> ~ 2 X D - *)K 1 = 1- 2 (st ) 2 —^F( 3. *) 


cos 


i2__ 


2.4.6.8 


(D - 3)(D ~ S)(D - 7 )h 1 = 1-3.5-JJff^ 2 (st) 3 -L^F( i, 4- *) 

'4' ’ cos 2 — 

2 


cos 


,1 


(D + 4 )(D - 3 )(D - 5)(D - 7)\ ! = 1 . 3 . 5 - “TTTTi 2 < st ) 8 :F ( ~ b ~b 4- *) 


COS 0 


e 


(D + 4)(D + 6)(D - 3).. (D 7)^s.i — I, 5 2.4 .... 12 2 (st ) 3 ^ F(--|, — 4 >:B ) 


(D - 4 )(D - 6)(D - 8)(D - 10)^ , = I- 3 -S -7 2 ( st ) 4 ^ F ( ~ i i» S - *) 

’ 4 * * * * cos 2 — 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjoumals.org/ at University of North Dakota on May 25, 2015 























1909MNRAS..69..6331 


June 1909. Development of the Perturbative Function. 645 


(D - 5)(D - 7).... (D - i 3 )^.i =I - 3 - 5 - 7 - 9 ^'g'y ^ 2 ( st ) 5 


F(-|, 6. x) 


COS 


2 _— 


(D - 6)(D - 8)-(D - i 6)6 6i j= 1.3.5 ... 1 i - 3 '| —~- 2 - 3 2 ( st ) 6 — ~ ¥> £> 7 - *) 

2.4*D • • • 24 


COS 




etc. 


etc. 


etc. 


We also have by direct differentiation of b 0t a 

0 


0 

D b n , = tan 0 tan— 

o-i o 


tan 2 - 


b 0 ’ 1 + cos 0 2 - x '> 


The formula for b Q1 is only inserted in the above table for sake 
of symmetry. It can be further transformed, as is well known, 
but the process used for its further transformation will not apply 
to the other formulae. 

Special formulae such as (D + 2)(D - i) b lt j = (D + 3)(D - 2) b Si 2 
also exist. 

If we apply some of these formulae to the case of Titan- 
Hyperion we have 

0 = 55 °’ 5949 22 

0 0 

log tan 2 0 = 0*32881604. log tan—- = 9*72196186. log cos 2 — = 9*89349547 

L ^ ^ 2 2 


(D - i)6 l i = - 5 -sin 0 tan 2 0 — J —- (i-^-x- ?-x 
’ 6 v c 


cos z 


8 64 


2 _ 5 ^,3. 


log series =9-99574456 
„ factors = 0*22686915 

(D - i)6 l 1 = 0*22261371 = 1*6696049 
The previous calculation gives 1*6696046 

Again, 

(D^i 5 D^ 7 iD-.os)A, 1 = 3 ^- S ( sin 6 ^.) 3 ( + Jj. s» - - 1 - 3 * - . . . 


5 12 


cos 


2___ 


16 128 


1024 


log series = 9*9895742 
,, factors = 0*6728478 

1*6624220 = 45*96434 
The tables give 45*96443 

With these coefficients we could proceed to form a theory of the 
motion of Hyperion as disturbed by Titan, but it is at least 
doubtful if the algebraical development of the perturbing function 
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646 Mr. Innes , Certain Coefficients in the Algebraical lxix. 8, 

converges in this case. K. F. Sandman, in considering movement 
in a plane (which is very nearly the case with Titan-Hyperion), 
shows that the development according to powers of the eccentricities 
is certainly divergent if 

a(i + 26 + e 2 + etc.)>(i — 2ej + ef - §ef + § e\ — etc.) 

(Ueber die Storungen .... Helsingfors, 1901). 

In the present case we have 

a 0*825 
e = 0*029 
e 1 = 0*104 

so that the inequality is 

o*874>o*8oi6 

and the development is divergent. In this case e 1 is not really an 
eccentricity in the planetary sense of the word; it is a perturbation 
due to the action of Titan; but this probably does not invalidate 
Sundman’s criterion. 

Newcomb and H. Struve have computed several of the coeffi¬ 
cients (.Astf. Papers , iii., and Poulkovo , xi 2 . pp. 267 and 283). 
Using H. Struve’s value of a (0*8250634) and Newcomb’s clear 
notation we have 

k x — ejPJ; \ 4 - efPJ; 1 + efPJ; f + etc. 

?Vi = i (7 + V)b Si i = +3*274193. 

Po;i=tV( ~ 358 ~ 37^ + ioD 2 + d 3 )& 3i i = + 3*534766. 
po .5 = ¥ i T ( I 6o72 + 286D - 942D 2 + i5l> 3 + 17I) 4 + D 6 )& 3( ! 

= + 262*4885 

k 2 = e\ PJl + ^Pfti + etc. 

P°’J = £(202 + 2 9 D + D 2 )Z> 6il = + 14*78513. 

PS; 2 - irV( - 415 7 2 “ 5360D + 113D 2 + 34D 3 + D 4 )fr 6i ! 

= - 112*1. 


The value of PS;2 is derived from Newcomb’s formulae; those 
of Le Yerrier give — 5062 as the coefficient of D& 6>1 . 

^ 3 = e?PS;3 + etc. 

Po;3 = t&(9 8 46+ 1415U + 66D 2 + D 3 )& 9(1 = +82*31306. 


P 0 = J(D + D 2 )fr 0il = +9*069466. 
p =|( 7 + U)(- 3 6 + D + D 2 ) 5 3( i = J( - 63 - 2D 

+ 8D 2 + D s )6 3 ,i = 


i 39 * 53 6 3 - 


N.B. *(&«_!, 3 + b i+h 3 ) = (D 2 + D - i*)b it ,. 

P 2 = £( 2 02 + 29D + D 2 )( - 36 + D + D 2 )& 6i L = 

- 7272 - 842D + 195D 2 + 30D 8 + D 4 )6 6(1 = + 1106^424. 
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’■ It will be noticed that the same coefficients repeat themselves 
in the last two quantities, but they are now operated on by 
(D2 + D-36). 

Lastly, the values of the coefficients may be required for a 
slightly different value of a, say a + Aa. If so, the figures given for 


a will be useful. Let us represent any coefficient D m ^ 
then by Taylor’s series 



fa + A«) =fa) + ff/(a) + i(^)*a s ffa) + etc. 


In this replace the operator 


a— by D which gives 
da 


a 2i!L = D(D-i) 

da* 

..^=1X1,-, KD-A 

etc. etc. 


and Taylor’s series becomes 

yia + Aa) =yip) + log ^1 + —^ vf(a) + | log ^1 + J D 2 J\a) + etc 


If-— is so small that its square is negligible we have 
a 

1 ( , Aa\ Aa 

l0e V + ^r~a- 

If the Newcomb coefficient D m b in be put in place of (a) we have 
putting l for log ^1 + —^ , 


For a + Aa For a 

TP%^ n = + + 

which enables us to pass readily from the coefficients calculated for 
a to those for a + Aa. 

Johannesburg, S.A. : 

2 6th March 1909. 


4 6 
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648 Mr. Robbins , Notes on Mr. Innes ’ Paper lxix. 8, 


Notes added by Frank Robbins. 

The two equations for b 0> 1 may be derived from equations (4) 
and (7) respectively as given by Yerhulst (loc. cit.) if it is 
remembered that his I) is the complete elliptic integral of the first 
class 

dcj> 

Ji - c 2 sin 2 <j> 


an.d his b is the complementary modulus Ji — e 2 . 

The second and more accurate formula may be slightly 
simplified—it easily becomes 



& o,i = 


y8(a| + a|)} 




but even so, it is somewhat troublesome, and the first formula would 
probably be adopted by most computers. It is certainly quicker 
even if a special value of N 1 is computed instead of being inter¬ 
polated from the table, page 649. 


Computation of K r 

The exact value of Nj is log 10 (1 + 2 q 4: + 2g 16 + 2q S6 + . . . .) 2 

or lo gio l 1 + 42 4 + 42 8 + ....)• 

If a is not greater than 0*85 it will be sufficient to stop at 4g 4 , 
and then it is easy to see from the equations 


a 2 + a? = 1 


2 


I - a\ 
2(1 -f- a}) 


l \ 9 


q = T + \T) + 1 \T) + ■ 

that N 1 = m log e (1 + 4g 4 + 4 q 8 + . . 
or Nj = 4 mqf very nearly 


•) 


Ni = 4 »»{ (ffj +8 (“i“) | m being 0-4342945. 


This is a very convenient working formula, and it can be com¬ 
puted for any given value of a in less than ten minutes. Zech’s 
addition table serves admirably here; the argument being the 
difference of the logarithms of the two terms within the brackets. 
This difference can be found mentally from the log of the first 

term ^ without the log 8 being known. It is merely 


necessary to add 0*9030900 to the mantissa of log 



and to 
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